Magnetic component of Yang-Mills plasma 
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Confinement in non-Abelian gauge theories is commonly ascribed to percolation of magnetic 
monopoles, or strings in the vacuum. At the deconfinement phase transition the condensed magnetic 
degrees of freedom are released into gluon plasma as thermal magnetic monopoles. We point out that 
within the percolation picture lattice simulations can be used to estimate the monopole content of 
the gluon plasma. We show that right above the critical temperature the monopole density remains 
a constant function of temperature, as for a liquid, and then grows, like for a gas. 
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It is well known, the properties of the Yang-Mills 
plasma turned out to be unexpected. In very brief, the 
plasma is similar rather to an ideal liquid than to a gluon 
gas interacting perturbatively Amusingly enough, 
many features of the plasma find their theoretical expla- 
nation in terms of a dual, or string formulation of Yang- 
Mills theories [2j which has been derived only in the limit 
of infinite number of colors and supersymmetry. 

It is a challenge to uncover dynamical picture behind 
the observations on plasma. In particular, equation of 
state has been studied in detail through lattice simula- 
tions both in case of SU (2) and SU (3) gauge theories Q. 
It turns out that global characteristics of plasma are not 
in contradiction either with perturbation theory or with 
the string picture, or with quasiparticle models and 
do not provide much insight into the plasma dynamics. 

Infrared-sensitive variables could be more helpful to 
identify specific degrees of freedom of the plasma. An 
example of such a variable is the viscosity which turns to 
be very low [f|. Other examples are the string tensions 
a and a of the spatial Wilson and 't Hooft loops, 

(Wspatiai) ~ exp(-ovl) , (if S patiai) ~ exp(-aA) , (1) 

respectively. Here A is the area of a minimal surface 
spanned on the corresponding contour. 

Thus, we come to consider light degrees of freedom 
of the plasma. On the theoretical side, our guiding ob- 
servation is that degrees of freedom condensed at T = 
form a light component of the thermal plasma at T > T c . 
Since the confinement of color in non-Abelian theories is 
due to magnetic degrees of freedom ((J @, S| the magnetic 
component is to be present in the plasma as well. 

In the string picture one postulates existence of elec- 
tric strings which can be open on external quarks (Wilson 
line) and of magnetic strings which can be open on ex- 
ternal monopoles ('t Hooft line). At T = 0, the electric 
strings have a non-zero tension, measured through the 
Wilson line. One of the earliest proposals for the mech- 
anisms of deconfinement is percolation of electric strings 
Generically, percolation means that there is large 



(potential) energy E and large entropy S which cancel 
each other. Thus, within this mechanism it is long strings 
which become tensionless and percolate. Lattice data [101 ] 
provides evidence for large entropy S = — (dF/dT) v of 
the electric strings, since the average of the Wilson line 
(W) ~ exp(— F/T) can directly be related to free en- 
ergy F = E — TS . The potential energy continues to 
grow above the phase transition, T c < T < 2T r , and the 
entropy reaches values of the order of 10, Ref. [1 01 ] . 

On the lattice, the phenomenon of strong cancela- 
tion between energy and entropy was discovered first 
for magnetic strings At zero temperature mag- 

netic strings percolate through the vacuum and pro- 
vide disorder which causes confinement. Percolating 
strings are directly observable on the lattice and are 
known as center vortices Q. Non-Abelian action asso- 
ciated with the strings is ultraviolet-divergent, Springs ~ 
(Area) s trings/« 2 , where a is the lattice spacing. The total 
area of the strings, on the other hand, is a finite quantity 
in physical units, ((Area) strings) ~ Aqcd^4<2- If there 
were no cancelation between the string energy and the 
string entropy, the area of the strings would be in lat- 
tice units as well. Note that there is no area law for 
the 't Hooft loop at at T — and therefore the tension 
of the magnetic string vanishes. Within large- 7V C dual 
formulations this relation is satisfied explicitly 0]. 

Quantum mechanically, the lowest string mode is 
tachyonic. The tachyonic mode corresponds to the 
monopole condensation [8|. Since there is only a single 
tachyonic mode of the string, the string percolation can 
in fact be projected into percolation of monopoles [111 ]. 

A link between percolation and field theory is provided 
by the polymer formulation of field theory in Euclidean 
space [l2fl . One starts with classical action of a free par- 
ticle, S c \ — M ■ L, where M = M(a) is a mass parameter 
and L is the length of trajectory. By evaluating the Feyn- 
man path integral for the particle propagator one learns 
that the propagating (physical) mass is in fact 
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where the constants are known explicitly for a particular 
(lattice, for example) regularization. The tachyonic case 
corresponds, as usual, to «^p hys < 0. The tachyonic mode 
is manifested as an infinite, or percolating cluster. To en- 
sure "iphys ~ Aq CD the bare mass parameter M(a) and 
the entropy factor, C en tr/o in Eq. ^ are to be fine tuned. 
This, crucial condition of fine tuning between energy and 
entropy is satisfied for the lattice monopoles More- 
over, percolation of the monopole at short distances can 
be understood within the approximation (|2|), while at 
large distances the properties of the monopole currents 
correspond to a constrained system Observation- 
ally, the constraint is that the monopoles live on two- 
dimensional surfaces, or strings, for a review see 
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It is a generic field-theoretic phenomenon that particles 
which are virtual at T = 0, are becoming real at finite 
temperature and released into the thermal plasma. In 
our case, there exists a tachyonic monopole mode which 
is to disappear at the point of the phase transition, with 
magnetic strings emerging into the plasma [la ] . In order 
to realize a tempting opportunity to check this picture 
from the first principles offered by Euclidean lattice sim- 
ulations, one should be able to distinguish between real 
and virtual particles in the Euclidean formulation of the 
theory. The problem is that even real particles in the 
Euclidean space are off-mass-shell. 

Let us first address the problem of detection of real 
particles in lattice simulations in case of a free scalar 
field (for a related discussion see Refs. [3, E3])- Since 
the monopoles are observed as closed trajectories on the 
lattice, it is appropriate to utilize the polymer represen- 
tation At zero temperature a typical monopole en- 
semble consists of a large number of finite clusters of 
the monopole trajectories corresponding to virtual par- 
ticles and one infinite cluster associated with condensed 
(tachyonic) monopoles. As temperature increases the in- 
finite cluster disappears since the monopole condensate 
vanishes at the deconfinement temperature. We point 
out that the monopole condensate melts down into real 
thermal particles, which contribute to magnetic content 
of the thermal plasma. 

In the imaginary time formalism the finite tempera- 
ture T is imposed via compactification of the time di- 
rection, £4 into a circle of length 1/T, and the points 
x = (x, X4 + s/T), s G Z, are identified. In the language 
of trajectories the integer number s has meaning of the 
wrapping number. It is obvious that properties of ther- 
mal particles are encoded in the wrapped trajectories, 
s ^ 0, and the virtual particles are non-wrapped, s = 0. 

The propagator of a scalar particle is given by 



G(x -y)oc 



~Scl[Px 



where the sum is over all trajectories P x . y connecting 
points x and y. Evaluation of the propagator at a finite 
temperature T is a straightforward generalization of the 



T = case [12]. The propagator in momentum space, 

g.(p) = J d 3 xe-*' x ) G(x,t = s/T) (3) 
is related to the thermal distribution of the particles 

1 g wr (p) 
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g m = Y,G s , G vac = Go, (4) 



given by the Bose-Einstein formula fx = l/(e"p/ T — 1) 
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Equation Q demonstrates that wrapped trajectories 
in the Euclidean space correspond to real particles in 
Minkowski space. However, the normalization to the per- 
turbative T = propagator, C* vac = 4/(a 2 oj p ), is awk- 
ward since it depends explicitly on the lattice spacing. 
Also, we do not expect in fact that the magnetic fluc- 
tuations of the Yang-Mills plasma - as probed by the 
lattice monopoles - correspond to free particles. As it is 
emphasized above, the monopoles represent only a com- 
ponent of the whole plasma and their properties are to 
be constrained by the environment. 

Moreover, one can explicitly show that the average 
number of wrappings s in a time slice of volume V^d is 
most directly related to the density of real particles 



p(T) = n m . = (\s\)/V 3d . 
where, in the case of the free particles, 
d 3 P iVd.f. 
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and [i is the (positively-defined) chemical potential con- 
trolling the particle number. A similar relation is known 
for a gas of non-relativistic scalar particles 17[ and goe s 
back to the Feynman's theory of A-transition in 4 He 18 1. 

The interaction of the monopoles with the environment 
is expressed, in particular, in the number of effective de- 
grees of freedom N^.i. and in the non-zero chemical po- 
tential fi to be discussed below. We will treat Eq. ([5]) 
as definition of the density of thermal particles in the 
imaginary time formalism. 

Thus, the density of the thermal particles p(T) corre- 
sponds to the vacuum expectation of the number of the 
wrapped loops. By measuring the number of wrappings 
(in the Euclidean space) one can learn density of real 
particles at finite temperature (in Minkowski space). 

A change of the character of the monopole trajecto- 
ries in Euclidean lattice simulations near the point of 
the phase transition has been observed and discussed in 
many papers, see, in particular, 0, [lj|. The change, 
indeed, is disappearance of the percolating cluster and 
emergence of wrapped trajectories. Quantitatively, only 
a part of the percolating (tachyonic) cluster goes into 
the wrapped monopole trajectories, while the other part 
is released into the vacuum as virtual particles. 
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We conclude that, qualitatively, there is little doubt 
that at T = T c there is transition of a part of the tachy- 
onic mode into degrees of freedom of the thermal plasma. 
This is a spectacular phenomenon by itself and is a proof 
of reality of a magnetic component of the gluon plasma. 

To be quantitative, we need a detailed lattice data on 
the wrapped trajectories. In case of pure SU{2) gauge 
theory there is relevant data in literature, see, in partic- 
ular, Refs. lj| The most detailed data refers, 
however, rather to the length density than to the number 
density ([5|) of the wrapped trajectories 
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where L wr is the total length of the wrapped trajecto- 
ries in the volume The expression for the density of 
the wrapping number which enters the relation ([5]) can 
be obtained from ([7]) by replacing the total length of the 
wrapped trajectories L wr by its projection on the time 
axis. The two expressions coincide in the limit of static 
trajectories. In reality, the approximation of static tra- 
jectories for the wrapped loops is reasonable. Also, it 
has been checked [19( that the density (O is independent 
on the lattice spacing and scales in the physical units, 
as is expected for the density of the wrapping number. 
Thus, we use the data [3| on the wrapped trajectories 
to estimate the density of real monopoles in plasma. 

The data indicates existence of two distinct regions 
in the deconfincmcnt phase: the first region covers the 
range of temperatures T c < T < 2T C , and the second 
region corresponds to the higher temperatures T > 2T C . 
In the first region the density of the thermal monopoles 
is almost insensitive to the temperature 0, [l9| 



p(T) m T 3 (T C <T< 2T C 
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The numerical value of the density - corresponding, ap- 
proximately, to 3.5 monopoles per cubic fermi - is com- 
parable to the density of monopoles in the monopole con- 
densate which amounts to approximately 7.5 monopoles 
per cubic fermi at zero temperature .20]. Thus, a siz- 
able fraction of the T = condensate is released at the 
transition region into the vacuum as thermal particles. 

It is instructive to compare the ratio (jSJ) with a limit 
of free relativistic particles. Then the ratio (JSj) at T = T c 
would be equal or smaller than a well-known number, 
C(3)/7r 2 w 0.12. Thus, the monopole medium just above 
the critical temperature is an order of magnitude denser 
than the ideal gas estimate. This fact along with the 
observation that the monopole density is independent on 
temperature allows us to suggest that in the region T c < 
T < 2T C the monopoles form a dense (magnetic) liquid. 

There is another radical difference between the mag- 
netic monopole constituents of the Yang-Mills plasma 
and free particles. For free particles, the wrapped tra- 
jectories are not static at all. Moreover, the density of 



the wrapped trajectories diverges at small lattice spac- 
ings a as p^ e ~ T 2 /q. On the other hand, the avail- 
able lattice data [l9( provides no indication of such a 
divergence and this can be explained only by specific 
constraints, encoded, for example, in chemical potential 
for the monopole trajectories, or in the constraints that 
monopole trajectories belong to surfaces, see also [la ]. 

Starting from T « 2T C the density of particles associ- 
ated with the wrapped trajectories grows, 



p(T) « (0.25 T) 3 (T>2T C ) 



(9) 



where we drop subleading terms in the asymptotic limit 
T — > oo. If we again compare ((9]) with the ideal gas case, 
then in the asymptotic limit the monopoles take much 
less than one degree of freedom 21]. 

According to the dimensional reduction arguments, 
valid at high temperature [IH non-perturbative physics 
is described by 3d magnetodynamics which corresponds 
to zero Matsubara frequency of the original 4d theory. 
In terms of the monopole trajectories the restriction to 
the zero Matsubara frequency implies that the wrapped 
trajectories become static. And indeed trajectories of 4d 
monopoles become static at higher temperature. Accord- 
ing to Refs. 23|, ]24| description of monopoles in terms 
of 3d, or static trajectories and 4d wrapped trajectories 
match each other at T sa 2.4T C . 

Matching with the dimensional reduction allows for an 
important consistency check of Eq. Namely, as far 
as the monopole trajectories are static and, consequently, 
non-intersecting they could be treated within a 3d the- 
ory [25| utilizing methods of (26|. Then the density ([5]) 
is equivalent to counting of monopoles trajectories. This 
procedure is true also in case of strong interaction be- 
tween the monopoles provided they are static. The di- 
mensional reduction implies that the monopole density is 



p(T) = C p g e 3d (T)K 
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logT/Aqcu) 



T>T C , 



(10) 

where C p is a temperature-independent parameter. Nu- 
merically, Eqs. ([9]) and (fT0"|) are compatible with each 
other within accuracy of the available lattice data. 

The temperature dependence exhibited by Eq pTJj) can 
be reproduced by Eq. ([5]) provided that there exists 
temperature-dependent chemical potential [2~fl ] 



p ~ Tlog.g 4 - 6 (T) ~ 3TloglogT/A, 



(11) 



which suppresses the monopole density (|6|) by the loga- 
rithmic factor, exp{-/i/T} - fff d (T) ~ l/log 3 (T/A). 

Thus, the evolution of the magnetic component of the 
Yang-Mills vacuum can schematically be represented as 



condensate 
(T < T c ) 



magnetic liquid 
(T C <T< 2T C ) 



gas 
(T > 2T C ) 



(12) 
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In the confinement region the magnetic monopoles con- 
stitute a colorless tachyonic state known as the monopole 
condensate while in the deconfincmcnt state they form at 
first a magnetized liquid and then a gas. In terms of the 
strings, percolation of both magnetic and electric strings 
at T > T cr ensures area law for both spatial Wilson and 
't Hooft loops, see © [27[. 

Because of the limited accuracy of the lattice data exis- 
tence of the chain (fT2|) is established rather on qualitative 
level. Further numerical studies seem to be well justified. 

In conclusion, let us emphasize analogy between phe- 
nomena in Yang-Mills theory with physics of superfluid- 
ity. In case of liquid helium there exist a superfluid and 
ordinary components of liquid. With increasing temper- 
ature particles from the superfluid component are trans- 
ferred to the ordinary- liquid component [28| . This is an 
analogy to the deconfincmcnt phase transition with vac- 
uum condensate vanishing and magnetic degrees of free- 
dom being released into the plasma around T « T c . In 
terms of the monopole trajectories the transition is from 
percolation in all directions to time-oriented trajectories. 
The total density of percolating and wrapped trajecto- 
ries remains approximately the same, in analogy with a 
superfluid. The "ordinary-magnetic-liquid" component 
might be responsible for low viscosity of the plasma. 

At higher temperature, T > 2T C the density of the 
magnetic component grows and approaches perturbative 
regime, p(T) ~ T 3 . This is analog of evaporation of the 
liquid. One might expect that beginning with tempera- 
tures T 2 T c the properties of plasma change gradually 
towards predictions of perturbation theory 
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